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FREE BOUNDARY PROBLEMS FOR TUMOR GROWTH: A VISCOSITY 

SOLUTIONS APPROACH 

INWON C. KIM, BENOiT PERTHAME, AND PANAGIOTIS E. SOUGANIDIS 


Abstract. The mathematical modeling of tumor growth leads to singular “stiff pressure law” limits 
for porous medium equations with a source term. Such asymptotic problems give rise to free bound¬ 
aries, which, in the absence of active motion, are generalized Hele-Shaw flows. In this note we use 
viscosity solutions methods to study limits for porous medium-type equations with active motion. 
We prove the uniform convergence of the density under fairly general assumptions on the initial 
data, thus improving existing results. We also obtain some additional information/regularity about 
the propagating interfaces, which, in view of the discontinuities, can nucleate and, thus, change 
topological type. The main tool is the construction of local, smooth, radial solutions which serve as 
barriers for the existence and uniqueness results as well as to quantify the speed of propagation of 
the free boundary propagation. 
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1. Introduction 


Motivated by models of tumor growth (see for instance the survey papers by Friedman [F], Lowengrub 
et al [LFJ“^]) and extending Perthame, Quiros and Vazquez [PQV], in a recent paper Perthame, 
Quiros, Tang and Vauchelet [PQTV] studied the limiting behavior, as m —)■ oo, of the solution 
(density) pm of the porous medium diffusion equation (pme for short) 

Pm,t-PmG{pm) in Qt X (0, T), (1.1) 

where il is either a bounded domain in K" or = K", with boundary and initial conditions 


Pm = Pl in dpQr ■= dV, x [0, T) and = Po.m on x {0}, 


( 1 . 2 ) 


satisfying 

0 < Pl < 1 and 0 < po < 1 if 11 is bounded or po G L^(IR") if H = K”. (1.3) 

Here 

Pm ■■= ^Pm-^ (1-4) 

is the pressure, v > 0, and G : K —>■ M is a smooth function, which describes the cell multiplication, 
satisfying 

G{pm) = 0 for some pm > 0 and G' < 0. (1.5) 

Using the pressure variable, (1.1) can be rewritten as 


Pm,t ^^'^{PmDpm ^Up^) — PmGij^m')^ 


a form which represents better the mechanical interpretation of the model with Vm ■= —Dpm the 
tissue bulk velocity according to Darcy’s law. 
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Note that, if, as m —?> oo, the Pm’s and pm’s converge respectively to p and p, then p will be nonzero 
only where p = 1. This indicates, that, in the limit, a phase transition may take place with an evolving 
free boundary between the tumor region (the support of p) and the pre-tumor zone (the support of 
1 — p). The convergence, as m —)■ oo, of the Pm’s and Pm’s has already been investigated using a 
distributional solution approach in [PQTV]. 

Here we study the asymptotic behavior of Pm and pm in the limit m ^ oo using viscosity solutions. 
This yields a different description of the limit problem, allows for more general initial data and yields 
pointwise information about the free boundary evolution, uniform convergence results as well as some 
quantified statements about the speed of propagation of the tumor zone. 

In order to state the result it is necessary to introduce the limit problem which we derive next formally 
following [PQTV]. We use the auxiliary variable 

Um •— Pra T Pm)-i (^’d) 


which is close to —pm + p(1 — Pm) for large m, and, recalling that (1.1) can also be written as 

m - 2 jDpmP 


Pm,t - (to - l)pmAp - \Dpm\ “ ^^Pm = (to - l)PmG(Pm) “ V 


m — 1 


(1.7) 


we find that Um satisfies 

[bm{Um)\t - = -VPmG(Pm) with b'^{u^) = --•. (1.8) 

mpm + V 

Assume next that, as m ^ oo, the Pm’s, Pm’s and Um’s converge respectively to p, p and u. Then, 
formally, we find 


p = u , v{\ — p) = where := max(M, 0) and u := — min(u, 0). (1-9) 


Letting to —)• oo in (1.7) and noting that {p > 0} = {p = 1} yields that p and p solve respectively 


- Ap{-,t) = G{p){-,t) in n{t) := {p{-,t) > 0}, 


( 1 . 10 ) 


and 

pt - vAp = pG(0) in {p < 1}; 

this last equation which can be rewritten in terms of p := iz(l — p) as 

Pt — vAp =—vpGiO) = [p — v)G{{)) in K" — r2(t). (1-11) 

Combining (1.9), (1.10) and (1.11) the limiting problem can be recast with h{u) := u+ as the elliptic- 
parabolic equation 

b{u)t — vAu = {b{u) — v)G{u~) in Qt- (1-12) 

In view of the Lipschitz continuity of b, (1.12) yields the free boundary condition 

drjU~^ = drjU~ on dn{t), (1-13) 

where t] is the (outward) normal at {x,t) € dfl(t), which provides an implicit motion law for the free 
boundary d{u > 0}. 

This problem but without the growth term has been studied by Alt-Luckhaus [AL] in the weak 
(duality) setting, by Carrillo [C] in the weak (entropy) setting and also by Kim-Pozar [KP] using 
viscosity solutions. The right hand side, however, plays an important role because it can generate 
pressure nucleation, thus generating a change of topology in the free boundary (in a zone where p = 0, 
p can grow and reach p = 1 thus generating a new island where p > 0). Viscosity solutions for porous 
medium equations were introduced by Caffarelli and Vazquez [CV], see also the book by Vazquez [V], 
and later by Kim [K] for the Hele-Shaw problem, for the Stefan problem and additional references the 
reader can refer to the book by Caffarelli and Salsa [CSj. 

Next we introduce the precise assumptions on the initial data Pm.o- We assume that 
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and, as TO —)■ c», 


there exists Mq > 0 such that := 


TO 




^ Mo 


i Pmfi —t Po uniformly, where po : “t [0,1] is continuous and 
Oq := {po = 1} is bounded domain with locally Lipschitz boundary, 
and, in addition, if = K" that 

Po G L\W^). 


The main result of the paper is: 


(1.14) 


(1.15) 

(1.16) 


Theorem 1.1. Let pm solve (1.1) with data satisfying (1.3), (1-14), (1.15), and (1.16), if ft = M", 
and define uo(x) := —wxfio + ~ PeXxo.^, where w be the unique solution of —Aw = G{w) in 

flo o,nd w = 0 on dflQ. Then, for all T > 0 and as m ^ oo, the Pm’s converge to 1 — v~^h(u) 
uniformly in Qt, where u is the unique viscosity solution of (1.12), with the initial and boundary data 
uq and v{l — pl) respectively. 


We remark that the Pm’s converge uniformly to U- as long as U- is continuous. It turns out, however, 
that U- may be discontinuous in time. This is due to the fact that (the positive part) b{u) solves a 
parabolic equation with a sink term, which means that h{u) can decrease to zero in the interior of its 
positive phase and nucleate a negative phase. Once the negative phase is created, the elliptic equation 
that needs to be satisfied in the negative phase leads to the jump discontinuity of u_ over time. We 
refer to Section 5 for a discussion about the propagation of the elliptic and parabolic phases. The 
Figure 1 illustrate the time discontinuities in pm, and Figure 2 shows and additional discontinuity in 
the Pm’s when viscosity is neglected (Hele-Shaw problem). 






Figure 1. Snapshots of the evolution of the density pm (blue solid line) and pressure 
Pm (black discontinuous line). The parameters are to = 20, v = .5. These figures 
illustrate how the pressure profile is building up when density approaches one. Be¬ 
tween the last two pictures, the density has continuously reached p = 1 in the center, 
while p has jumped discontinuously. 

Theorem 1.1 is proved by showing first that the Mm’s converge to u as to —> oo (see Corollary 4.4). The 
convergence results for the Pm’s then follow from the definition of Um and the continuity properties 
of h{u). 

As already mentioned above we use viscosity solutions which are based on appropriate choices of test 
functions. In the problem at hand these will be radial smooth solutions which are defined in Section 2. 
The main step in the proof of Theorem 1.1 is the following: 

Theorem 1.2. Assume that cf is a classical radial solution of (1.12). Then there exists a family of 
radial sub- and solutions super-solutions of (1.8) converging, as m ^ oo, to 4 >. 
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Figure 2. Snapshots of the evolution of the density pm (blue solid line) and pressure 
Prn (black discontinuous line). The parameters are m = 20, u = 0 (not treated in this 
paper). Compared to Figure 1, we observe also a discontinuous behaviour of the 
density. 

Organization of the paper. In Section 2 we introduce the notion of viscosity solution of the limit 
problem and discuss its existence and uniqueness which rely on the results of [AL] and [KP]. In 
Section 3 we prove Theorem 1.2. We used the results of Section 2 and Section 3 to prove Theorem 
1.1 in Section 4. In Section 5 we present estimates on the speed of propagation of the interface. In 
the Appendix we touch upon the proof of the comparison for viscosity solutions to the limit problem. 

Notation and Terminology. A nonempty set E C M" x K is called a parabolic neighborhood if 
E = U n (K” X (0, t]) for some open set U C M" x M and some r S (0, oo) and dpE := E \ E is the 
parabolic boundary of E. Whenever we refer to a parabolic neighborhood, we assume U and r are 
known. For an open subset A of K" x K and a continuous function u : A —)■ K, where A is the closure of 
A, {u > 0}, {v = 0} and {v < 0} are respectively the sets {y G A : v{y) > 0}, {y G A : v{y) = 0} and 
{y G A : v{y) < 0}. For a set A, we write XA for its characteristic function. If A, B are subsets of i?^, 
A + B := {a + b : a G A, b e B}. Given a set ?7 C USC{U) and LSC{U) are respectively the sets 
of upper and lower semicontinuous functions on [/; if we are dealing with bounded semicontinuous 
functions we write BUSC{U) and BLSC{U). If u : 17 —?> M is bounded, u*’^ and u* p denote 
respectively its upper and lower-semicontinuous envelope. We write Bp{x) for the open ball in K” 
centered x G K" with radius i? > 0. The duality map between a Banach space and its dual is denoted 
by < •,• >. If A is open subset of M” x M, Cl’l(A) denotes the space of functions on A which are 
continuously twice differentiable in x and continuously differentiable in t; similarly, if B is an open 
subset of K^', Cl.{B) denotes the space of 1-times continuously differentiable functions on B; when 
there is no confusion we omit the dependence on A and B. We say that a constant is dimensional if 
it depends only on the dimension. We denote by r) the gradient at (^, r) of ^ : 17 —>■ K, where 

U is an open subset of K" x K and, D4A{^,t) := lim (j, D(p{x,t)- For a family of functions 

fm'.rix (0, oo) -A K, liminf*/^ and lim sup*denote the lower and upper semi-continuous limits, 
that is liminf*/m(a;, 1) := limr_s.oinf{/m(T-l-y, t-l-s) : (x + yA + s) G Clx (0,oo), |y| + |s| <r,m> r~^}, 
and lim sup*/m(a;, 1) := limr_>o sup{/m(a; + y, 1 + s) : (a;-|-y,t-|-s) G x (0, oo), |y| -I- |s| < r, m > r~^}. 

2. Viscosity and weak solutions of the limit problem 
We introduce the notions of viscosity and weak (distributional) solutions to the limit problem 

{ b{u)t — v/S.u = {b{u) — u)G(u“) in Qt, 
u = uq on n X {0}, 
u = g on dil x [0, T], 


( 2 . 1 ) 
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and we discuss their existence and uniqueness. While we rely on the theory of weak solutions to obtain 
uniqueness results, viscosity solutions are used to prove pointwise convergence results in Section 4 as 
well as to derive information about the evolution of parabolic and elliptic phases in Section 5. 

Viscosity solutions: definition and comparison principle. Following [KP] we define viscosity 
solutions of (2.1). We begin by introducing the class of test functions which are classical sub-and 
super-solutions of the problem with the specific properties stated in the next definition. 

Definition 2.1 (Classical sub- and super-solutions). E be a parabolic neighborhood. Then (f € C{E) 
is a classical subsolution of (1.12) in E if 

(i) if G Clj{{ip > 0}) n Cl]l{{ip < 0}), 

(ii) {(^ = 0} C d {93 > 0} n d {</5 < 0} and \Dip^\ >0 on {(f = 0}, 

(Hi) b{if)t — v/S.ip < {b{ip) — iy)G{ip~^) on {f > 0} and {f < 0}, and 
(iv) \D(p'^\ > \Dip~\ on {if = 0}; 

if is a strict classical subsolution if the inequalities in (Hi) and (iv) are strict. Supersolutions are 
defined similarly by reversing the inequalities in (Hi) and (iv). 

Next we define the viscosity solutions to (1.12). We begin with the notions of (semicontinuous) sub-and 
super-solutions. 

Definition 2.2 (Viscosity sub- and super-solutions). We say u G BUSC{Qt) (resp. u G BUSC{Qt) 
is a viscosity subsolution (resp. supersolution) to (2.1) in Qt, if 

(i) u(-,0) < uo (respectively u(-,0) > uq) on Tlx {0}, u < g (resp. u> g) on dTl x [0,T], 

(ii) u < if (resp. u > if) on E for any strict classical subsolution (resp. supersolution) if on any 
parabolic neighborhood E C Qt for which u < if (resp. u> if) on dpE. 

We say u has initial data uq if u converges uniformly to uq as t ^ 0. 

Given that we do not expect to have continuous but rather lower semicontinuous solutions to (2.1) 
(recall the earlier discussion about possible nucleation), we define the notion of viscosity solution for 
such functions. 

Definition 2.3 (Viscosity solutions). A function u G BLSC{Qt) is a viscosity solution of ( 2 . 1 ) in 
Qt if its upper semi-continuous envelope is a viscosity subsolution in Qt and u is a viscosity 

supersolution in Qt- 

Next we introduce conditions on the data of (2.1) that are necessary in order to have a well-posed 
theory. 

For the lateral boundary we assume, for simplicity, that 

g is continuous and positive. ( 2 . 2 ) 

For the initial condition, if Tl is bounded, we assume that uq : t K satisfies 

Mo G <^(^2), Mo < 0 on dTl and — Amq = G(mo) in {mq > 0 }, (2.3) 

and 

F(mo) := i9{mo > 0} is locally a Lipschitz graph. (2.4) 

When T1 = K", we require, in addition, that 

{M(a:, 0) < 0} C Bp{x) for some R> 0 and m(-,0) — u G L^(M"). (2.5) 

Note that, in terms of p = m^ and p = v~^u- +1, the above conditions imply that the pressure phase 
is initially bounded and the density is integrable. 

We recall that we use only strict sub-and super-solutions as test functions. As a result it is possible 
to narrow the choice of the parabolic neighborhood E in Definition 2.2 and to use, instead, parabolic 
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cylinders of the form Q' = 0! x (^ 1 ,^ 2 ] C Qt, with 17' C n having a smooth boundary. Indeed, 
suppose that E C Qt is a parabolic neighborhood, (/j is a strict classical supersolution on E, u < (p 
on dpE and u > (p at some point in E. Define r := sup{s : u < (p on E D {t < s}} G K. The set 
A := {x : (x, r) € E,u > (pj is compact and, therefore, d := distance(A x {r}, dpE) > 0. Consider 
the parabolic cylinder Q' = {A + Bci/ 2 ) x (r —d, t]. Clearly Q' C E and u < (p on dpQ'; the boundary 
ot A-\- B 4/2 can be easily regularized. 

The proof of the comparison principle follows along the lines of the one in [KP] the only difference being 
the contribution of the source term {b{u) + i')G(u), which is not present in the equation considered 
in [KP]. The needed modifications are presented in the Appendix. Here we note, that, although, as 
discussed earlier, the source term in the negative (elliptic) phase of u may give rise to nucleations, 
this does not hinder the argument in [KP]. Indeed nucleations can not occur at the contact point of 
two “regularized” solutions which are initially strictly ordered (see Appendix). Hence we have the 
following theorem, which corresponds to Theorem 3.1 of [KP]. 

Theorem 2.1. Assume (2.2), (2.3), (2.4) and, if ft = K", (2.5) and let u G BUSC{Qt) and v G 
BLSC{Qt) be respectively a sub- and a super-solution of (1.12) in Qp for some T > 0 . If u < v on 
the parabolic boundary of Qt, then u < v in Qt- 

Note that in the above theorem we require that u and v are strictly separated on the parabolic 
boundary of Qt- For this reason. Theorem 2.1 yields, for given data, only the maximal and minimal 
viscosity solutions. Removing the condition of strict separation remains an interesting (and important) 
open question for most of free boundary problems. For the problem at hand, we are able to borrow 
the theory of regular weak solutions of [AL] to study the stability and uniqueness properties of the 
viscosity solutions. 

Weak solutions: existence and comparison principle. We recall the notion of regular weak 
solution to (2.1) and then state the theorems in [AL] which concern the existence and stability of such 
solutions. The difference in the proofs due to the source terms are discussed in the Appendix. We 
note that we do not try to recall the full generality of [AL] but we modify definitions and statements 
to apply to the specific problem we study here. Finally, we remark that [AL] only considers bounded 
domains D. 

We begin with the definitions of weak and regular solutions. The former is a notion of solutions based 
on duality (integration by parts) while the second, as the name indicates, requires more regularity (in 
time). 

Definition 2.4 (Weak solutions), u G L^(0,T;iL^(fl)) with b{u) G L°°{0,T;L^{n)) and dtb{u) G 
L'^{ 0 ,t; is a weak solution to to (1.12) in Qt if, for every test function C. € Lf(Q,T-,HQ{Vt)), 

f < dtb{u), ( > dt-\- f f Du-Dfdxdt= f ( {b{u) — p)G{u~)fdxdt. (2.6) 

Jo Jo Jn Jo Jq 

Weak sub- and super-solutions are defined with the corresponding inequalities replacing the equality 
above. 

Definition 2.5 (Regular weak solutions). A weak solution u is regular if dtb{u) G L‘^{Qt)- 
The following corresponds to Theorem 2.2 of [AL]. 

Theorem 2.2. Assume that the bounded domain 17 has Lipschitz boundary. Ifu and v are respectively 
a weak subsolution and supersolution of (1.12) in Qt and, in addition, dt{b{u) — b{v)) G L^{Qt) and 
u < V on the parabolic boundary of Qt, then u < v in Qt- 

The proof parallels the one of Theorem 2.2 in [AL]. It is rather simple, and, for the reader’s convenience, 
we present it next. 
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Proof. Fix (5 > 0 and let tp^z) := min[l,max(0/(5,0)]. Applying the definition to u and v with 
C, := 'tp{u — v) and using the fact that G is decreasing and u_ < if u < m, we find, for some C > 0, 

fo In 9t{b{u) - b{v))tp{v - u)dxdt + f fg X{o<u-v<s}ld^(u - v)l^dxdt 

^ fo fnlX{o<u-v<s}(b(u) - u)G(u_) - (b(v) - iy)G(v_)]dxdt 

= Jo fnlX{o<u-v<s}(b(u) - b(v))G(0) + u(-G(m_) + G(v_))]dxdt 
^ /o- b(v))+G(0)dxdt. 

Next we let 5 —>■ 0. The parabolic term on the left converges to f^(b(u) — b(v))+(-,t)dx, and, thus, 
Gronwall’s lemma yields that b(u) < b(v). 

Using this last information in the previous inequality, we find that, for any d > 0, b(u — u) = 0 in 
{0 < u — u < d}, which, in view of the fact that u < v on the parabolic boundary of Qt, yields the 
claimed comparison. □ 

The next results assert the existence and uniqueness of regular weak solutions. 

Theorem 2.3. Assume that the bounded domain has Lipschitz boundary. Then there exists a 
unique regular weak solution u in Qt for a given boundary data ud G H^{fl)). Moreover, 

b{u) is continuous. 

The uniqueness follows from Theorem 2.2 while the continuity of b{u) is a consequence of the main 
result of Di Benedetto and Gariepy [DG]. The existence is based on the Galerkin approximation and 
parallels that of Theorem 2.3 in [AL] where we refer for the details. 

Weak solutions to (1-12) are more flexible than viscosity solutions with respect to the range of the 
boundary data. On the other hand, we do not see how to prove directly the pointwise convergence of 
solutions of (1.8) to the solutions of (1-12) without going through the viscosity solutions framework. 
Indeed the convergence proof depends on the facts that we allow semi-continuity for viscosity sub- 
and super-solutions, and that, even with such low regularity, there is a comparison principle. 


Uniqueness and existence of viscosity solutions for bounded domains. We discuss now the 
uniqueness for the viscosity solutions and we argue separately depending on whether Xl is bounded or 
not.. 


Bounded domain. The first step to establish the uniqueness is that viscosity solutions evolve con¬ 
tinuously from the initial data. This is quantified in terms of the growth in time of the distance 
between the parabolic phase and dft \ flo, the initial position, where dflo ■= {u-o < 0}. 

Lemma 2.4. Assume (1.5), (2.2), (2.3) and (2.4) and let u be a viscosity solution of (2.1). There 
exists some sufficiently small to such that, for all 0 < t <to, 

d{x, d{uo > 0}) < t^^^ for any x G d{u{-,t) > 0}. 

Proof. To prove the claim we construct suitable super- and sub-solution barriers. We begin with the 
former. 

Fix £ > 0 and let S(t) be the 


S(t) := {a; : d(x, \ flo) < £ + 
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Let w'^ and w~ be respectively the solutions to 

= Ar(;+ in {(x,t) : x € S(t), 0 < t < 1}, 

= pl on dfl X [0,1], 

< 

= 0 on {(x,t) : x G dS(t), 0 < t < 1}, 
w+(-,0) = (uo)+ on S(0), 

and to the elliptic equation in (1.12) set in the complement of S(t) with zero boundary data, and 
define w := w~^ — w~. 

Since S(t) has the exterior ball condition, for small t > 0, we have 

\Dw-\{-,t) > e + on i9I](t). 

On the other hand, it follows, using as barriers, if necessary, shifted versions of heat kernel as barriers, 
that the restriction of \Dw+\ on 9S(t) converges exponentially fast to zero as t —)■ 0. 

Hence, we can choose a sufficiently small to > 0 such that, for all 0 < f < toj 

\Dw+\{-,t) < \Dw-\{-,t) on dT,{t). 

It follows that w is as a supersolution for (1.12). This yields that d{u{-,t) > 0} cannot expand faster 
than order of 

Similarly one can construct a subsolution barrier, based on the fast-decreasing set 

E(t) := {x : d{x,^lo) > 

Since the arguments are similar we omit them. □ 

Using the weak theory described above, we can prove, using arguments similar to the ones in Section 5.1 
of [KP], the following theorem. 

Theorem 2.5. Assume (1.5), (2.2), (2.3) and (2.4). There exists a unique viscosity solution u to 
(2.1). Moreover, u = (u*)*, and u coincide a.e. with the regular weak solution u of (1.12) in Qt- In 
particular, b{u) is continuous. 

Unbounded domain H = M". It can be checked easily, using a suitable radial barrier, that the 
negative phase of any viscosity solution of (1.12) with initial data uq satisfying (2.3), (2.4) and (2.5) 
is contained, for 0 < t < T, in with Ri = Ri{T). 

Let Un be the viscosity solution of (2.1) in Hflj+„(0) x [0,T] with initial data uoX{|a;|<_Ri+n} and 
boundary data —p — E on di?fl,^+„(0) x [0,T] respectively. It is immediate that the u„’s are mono- 
tonically increasing with respect to n, and, thus, converge pointwise to some w : M” x [0, T] —)■ M. 
Furthermore, the it„’s solve a uniformly parabolic equation outside of and, thus, are smooth in 
{i?i < |a;| < i?i -I- n} X (0,T]. We may, therefore, conclude that the convergence of the u„’s to w is 
locally uniform in {|a:| > i?i -|- 1} x [0, T) and, moreover, that w is smooth , since it solves a parabolic 
equation in {|a:| > -I- 1} x [0,T). 

Now we can use standard stability arguments about viscosity solutions to show that w is the unique 
viscosity solution of (2.1) in {|a;| < i?i -b 1} x [0,T] with the initial and boundary data uq and w on 
{|a;| = -b 1} X [0,r] respectively. It follows that w is a viscosity solution of (1.12) in M" x (0,T]. 

Theorem 2.6. Assume that uq satisfies (1.5), (2.3), (2.4) and (2.5). Then (1.12) has a unique 
viscosity solution in K" x (0,T] with initial datum uq. 

Proof. The barrier argument described above yields that any viscosity solution u to (1.12) in M” x 
(0, T], with initial condition satisfying the assumptions in the statement and has, in addition, compact 
non-positive phase, stays positive outside of a compact set, and, thus, solves the parabolic equation 
in (1.12) in {|a;| > R} x (0,T] for some i? > 0. 
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In particular, it follows that u and b{u) are smooth when |a;| is large, and b{u) — vG L°°(0, T; L^(IR")). 
Moreover, u is a viscosity solution of (1-12) in {|a;| < + 1} x (0,T] with smooth lateral boundary 

data. It follows that u is also the unique regular weak solution of (1-12) given by Theorem 2.3. 

This last claim and the decay of u for large |a;| imply that dth{u) G x {0,T)). The proof of 

Theorem 2.2 yields that u is the unique viscosity solution of (1.12) with initial data uq. □ 

3. The construction of radial barriers and the proof of Theorem 1.2 

We prove here Theorem 1.2. Let (phea radial classical subsolution to (1.12), as defined in Definition 2.1 
in the domain {ri < |a;| < r 2 } x [0,r], where 0 < ri < r 2 < oo, and, to simplify, we assume that the 
elliptic phase is contained in . 

We recall that this means that (p is smooth in its positive and negative phase and, for some a G 

Ci([ 0 ,r];K), 

< G{(j)~) in {()) < 0} = {{x,t) : ri < |a;| < a{t)}, 

4>t - vAp >{4>- u)G(O) in {()) > 0} = {(a:, t) : a{t) < |a;| < r 2 }, 

0 < |D^+| < \D4)~\ on {{x,t) : \x\ = a{t)}. 

We will perturb p to construct subsolutions to (1.8). It is, however, important to remark that these 
subsolutions do not have to be, actually will not be, smooth across their interface. Discontinous subso¬ 
lutions to the equation (1.1) can be defined in the viscosity sense in a fashion similar to Definition 2.2 

To keep things simple, we assume that (i) we have equality instead of inequality in the elliptic and 
parabolic equation in each phase of (p, and (ii) the domain is K" x [ 0 ,oo) and the positive (elliptic) 
phase of (p is {|a;| < a{t)}, that is we assume that there is no ri and r 2 = oo. A minor modification 
of the arguments presented below yields the general case. 

We introduce the notation 

Po := p- and po := 1 + v~^P+, 
and note that, for each t G (0,T], po{-,t) solves 

-Apo = G{po) in {|x| < a{t)}, 

Po = 0 on {ja;|=a(t)}, 

while Po satisfies 

pt — vAp = pG(0) in {{x,t) : \x\ > a{t)} 
with the free boundary condition 

\Dp~\ = iy\Dpo\ > \Dpo\ = lop'll on {{x,t) : \x\ = a{t)}. (3.3) 

Our first claim is: 

Lemma 3.1. There exist subsolutions Um to (1.8) which converge, as m ^ oo, uniformly to p, and, 
hence, there exist viscosity subsolutions pm of ( 1 . 1 ) such that, as m ^ oo, the Pm’s and the associated 
pressure variables Pm’s converge uniformly to po and po respectively. 

Proof. Consider the one-to one functions 

<i>(p) := -p^ + u(l - Pm) and ^{p) := $((!!i^p)i/™-i) ( 3 . 4 ) 

m 

and recall from the introduction that 

Um = $(Pm) = ^(Pm)- 

It is then immediate that, if Pm is bounded, then 


(3.1) 

(3.2) 


(3.5) 
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The heuristics in the introduction suggest as a possible way to prove the claim to perturb po and 
i/(po — 1 ) from the regions {po > 0 } = {|x| < a(t)} and {pQ < 1 } = {|x| > a(t)} and patch them 
together to construct a supersolution Um of ( 1 - 8 ), which then would imply that p := is a 

supersolution of ( 1 . 1 ). 

The barrier will be of the form 

-Um,i{\x\,t) in {|x| < a(t)}, 

Um, 2 (|a;|,t) in {|x| > a(t)}, 

where —Um,i and Um ,2 are classical subsolutions to ( 1 . 8 ) in their respective regions and at r = a{t) 
satisfy |dlui| < |Zlu 2 |. This inequality, which follows from (3.3), will prevent any smooth functions 
from crossing Um from below at {r = a(t)}. It follows that Um is a viscosity supersolution of (1.8) in 
the entire domain. 

We begin with Ui^m- Let Ura{-,t) solve 

-A'u„i(-,t) = G(pm) +/m(|x|) in {|x| < a(t)} and Urnia{t),t) = 0, 

where 

fmir) := Aoiy~^X{po<m-^/i}{r) + and Pm ■= (3.7) 

with Aq > 0 an independent of m sufficiently large constant to be determined below, and observe 
that, in view of the form of the equation above, Um is spatially radial and 

(^j f) — ’^m ( 7^ 5 f) ■ 
a[t) 

Note that, as m —>■ oo, / —>■ 0 in L^([0,a(t)]), since \Dpo\{a{t)) ^ 0. Using that p is radial, we find 
that, as TO —)■ c», u —)■ po in the G^-norm in {|x| < a{t)}. Since |71po| > 0 and u = po = 0 at |x| = a{t), 
it then follows that 

Po(x) = u{x) + o(|x — a(t)l) near |x| = a(t). 

In particular, in view of (3.6), we have 

fir) > ^oJ^”A{ii<|m-i/ 3 }(?’) + . (3.8) 

Next we define 

Ump := Mm - Cm with Cm '■= ^ {v/m^), 
and remark that, for any fc > I, 

— InTO < Cm < ^. (3.9) 

m 771^ 

The aim is to show that, for Um,i, Pm,i ■= ^~^i—Um,i) and Pm,i := ^~^(~Mm,i), 

(Mm,i)i - impm,i + h')Aum.i > {mpm,! + J^)piG(pm.i) m {(x,t) : |x| < a(t)}, (3.10) 

which yields that Um = —Um,i is a subsolution of ( 1 . 8 ) in the region {{x,t) : |x| < a{t)}. 

To show (3.10), note that, since Um,iix,t) = Um,i (aP)’ 1 ), it follows that, for sufficiently large to, 
|(Mm,i)t| = |(Mm)t| < |(x/a(t))'ZlUm| < G|Zlpo| = 0(1) in {(x,t) : \x\ < a{t)}. 

Also note that, in view of (3.9), the difference pm —Pm,i is of order of Cm if Pm > _ jn particular 

ifpm.i < then Pm < and, therefore, in view of (3.8), we find 


-{mpm,l + l2)AUmA = -{mpm,! + >^)AUm 

> {rnpm,! + J^)[G(pm) + + (mpm,! + ’^)ifir) - m~'^A) 

> [mpm,! + v)G{pm) + + (TOpm,i + v)ifir) - to“1/3), 

> {mpm.l + J^)Pm,T~^0(pm,l) “ iUm,l)t, 
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provided that m and ars sufScienty large. In particular the second inequality follows from the fact 
that, in view of (3.8), vfmiPm) > if Pm < while, otherwise, Pm,i > thus, 

TnPm.i > —)■ oo as m —>■ oo. 

It follows that (3.10) holds in {{x,t) : \x\ < a(t)}. 

We next define Um .,2 ■= ^(pm), where pm is a perturbation of po solving 

Pm,t - ’^‘^Pm = PmG{0) - j] — {{x,t) : |ai| > a{t),t > 0} 

with 

/Om(-,0) = Po(-,0) - 1 +$“^(cm) and $(p^) = on {|a:| = a(t)}. 

Note that 

^~^{cm) ^ (— lnm)^'^'"“^ —)■ 1“ as m —>■ oo, 
m 

therefore, 

Pm = (^)^/'""^ ~ 1 on {|a:| = a(<)}. 

Thus the Pm’s converge to po in the norm in E and in C], up to the boundary. Hence, since 

|Z?po|(a(t), t) > 0 and po is strictly away from 1 in S, for sufficiently large m we have p™“^ < vjw?, 
and, thus, 

A(p™) = mp™“^Ap„ + m{m - l)p'^~'^\Dpmf ^Ap„ + —\Dpmf = (3.11) 

m 

Then, for sufficiently large m, 

Pjn,t - A(p™) - i^Apm < Pm,t “ uAp^ + 0{m~^) < p^G(O) in S, 
and, therefore, Um ,2 is a supersolution to (1.8) in E. 

To conclude we check that the spatial gradients of Um,i and Um ,2 are ordered in the right order at the 
patching location |a;| = a{t). This follows since, in view of (3.3) and the G^-convergence of Um to po 
in {|a;| < a(t)}, 

\Du,n.l\ = \DUrn\ = (1 + o{l))\Dpo\ < (1 - o(l))u|£)po| < \DUjn, 2 \ On {\x\ = a{t)}, 
where the second inequality is due to (3.3) and the last inequality follows from the fact that 

|Zlp™| < m-2|i^p„| when p™"! = u/m^. 

□ 


Next we choose a radially monotone classical supersolution (j) of (1.12) and define po and po as before. 
Note that in this case we have 

|I1m+| = u|Ilpo| < |Tipo| = \Du~\ on {{x,t) : |a;| = a{t)}. (3-12) 

Lemma 3.2. There exists a supersolution Pm of (1.1) such that, as m ^ oo, Pm and Pm converge 
uniformly to po and po respectively. 


Proof. The argument parallels that of the proof of the previous lemma. 
Um{x,t) = Um{\x\,t) aS 


Um{r,t) 


-Um,i{r,t) in {r < a(t)}, 
'^ 771 , 2(^7 t), in {r > aft)}, 


with —Um,i and Um .2 solving, in the respective region. 


Um,t - {rnp^ + v)Aum > -p{mp^ + ty)G{pm), 


We define the barrier 


(3.13) 


and 


\DUm,l\ > \DUm,2\ Oli r = a{t). 
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This makes Um a (nonsmooth) supersolution of (1.8). 

We define Um,i = Um — Cm, where Cm = is as before and Um solves 

-Aum = G{pm) - fmix) in {|a;| < a(t)}, with Um,i(a(t)) = 0, 
with fm as in (3.7) in {|ai| < a(t)}. 

It follows that converges to Po{-,t) in the norm in {|x| < a{t)}. Moreover, since = 

0 (1), a straightforward computation, similar to the supersolution case, yields that 

Um,l,t - {mpm.l + l^)AUm,l < PmA{mpmp + l')G{pm,l), 
and thus (3.13) is satisfied. 

Next let Um ,2 '■= ^(Pm), where pm is a perturbation of po solving 

Pm,t - vApm = /OmG(O) + jn E := {{x,t) : |a;| > a{t),t > 0}, 

with Pm{-,0) = Po{-,0) - 1 + and ^{pm) = Cm on {|x| = a(t)}. 

Observe that, as m ^ oo, Pm converges to po ia the —norm in E and in the G,),—norm up to the 
boundary, and, since |I?po|(a(t), t) > 0 and po < 1 in E, for sufficiently large m, we find 

pZ~^ < Pm"^(a(t)) < = v/m^ in E. 

Thus it follows that 

HpZ) = mpZ~^^P-m + m{m - l)pZ~‘^\Dp^\‘^ - ^Ap^ + —|Op„p = 0 (to"^), (3.14) 

m 

and, for sufficiently large m, Um ,2 is a supersolution of (1.8), in view of the fact, that in this case 
Pm,t - A(p™) - vApm > Pm,t “ vAp^ “ 0(to"^) > pG(0) in E. 

Lastly, in view of (3.12), we have on {{x,t) : \x\ = a(<)} 

\Dum,i\{a{t),t) > (1 - o{l))\Dpo\{a{t),t) > (1 + o{l))iy\Dpo\{a{t),t) > \Dum, 2 \{a{t),t), 
which is the correct order for Um to be a subsolution to (1.8). 

□ 

Combining the two lemmas above yields the proof of Theorem 1.2. 

4. Convergence in the general case 

We use the results of the previous section to show the general convergence result, that is Theorem 1.1. 
We only consider the case fl = K”, since the arguments for bounded are similar. 

Let 

u := limsup*Um and v := liminf*Um. (4.1) 

We have: 

Theorem 4.1. u is a subsolution and v is a supersolution to (1.12). 

Proof. We will only show that u is a subsolution of (1.12). The proof of the other claim is similar. 

If u is not a subsolution, there must exist a classical strict supersolution ^ of (1.12) in a parabolic 
cylinder C := Br{xo) x such that u < (j) on the parabolic boundary of C and u crosses (j) from 

below for the first time at t = tQ. Let (xo,to) be where u — cf take its nonnegative maximum. 

Next we perturb (/) so that it becomes slightly smaller and thus u — 4> achieves a positive maximum in 
C. For instance, we may replace (j) by the inf-convolution 
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for sufficiently small r > 0 , and solve the corresponding elliptic and parabolic problem in the positive 
and negative phase of (p. Note that this perturbation preserves the supersolution condition < 

on the zero level set. 

The definition of u yields a subsequence, that we still denote by m, along which, for sufficiently large 
Ri) Um — (p has a positive maximum in C with Um < (f> on the parabolic boundary of C. Since Um is 
continuous, this means that Um = <^ at for the first time in C with ti < 

Let (yojSo) be a limit point of We claim that 

(yoj So) G > 0}- (4.2) 

To see this, suppose 4>{yo,so) > S > 0. Then the same holds at {xm,tm) for large enough m. It 
follows from the definition of Um that, as m —>■ oo, mpm —^ 0 uniformly. This and the fact that p 
is a supersolution of the parabolic equation in ( 1 . 12 ) yield a contradiction by a maximum principle 
argument, since, in view of ( 1 . 6 ), at {Xm,tm) 

4>t - < Um.i - V^Ura 

< mpm^Um - mPmPmG{pm)vpmG{Pm) 

< -V PmG{Q) + o{l) 

< (<()-u)G(0)+o(l). 

Note that in the third inequality of above computation we have used that 

rnpm^UmixmGm) < 'mpm^<P{^m,tTn) and, in view of the regularity of </> in the positive phase, 
limTO_>oo mpm^4>(xm, tjn) = 0. In the last inequality we have used that —vpm = Um — u. 

Similar arguments yield that the case {yo, sq) G d{(j) > 0}. 

For the rest of the proof we consider the subsequence of m’s along which the {xm^tm)'^ converges 
to (j/OjSo)- Since p has boundary, there exists a spatial ball Br{zo) C {^(•,so) < 0} which 

touches the free boundary at yo- 

Now we approximate p hy a radial supersolution tp of (1.12) in the space-time domain 

E := {a; : ri < |a; - zo\ < r2} x [sq - G sq], 

for some ri < r < r 2 - Roughly speaking the construction amounts to taking Taylor’s expansion of p 
at (j/ 0 )So) in each phase, up to the first order in space and time, and constructing a radial function 
with them. 

Let s be the outward normal velocity of the free boundary d{p > 0} at {yo, Sq), pick e > 0 small, let 
r{t) := r + a{t — Sq) for a = s + e so that Br(t){zo) C {p{-,t) < 0} for ti < t < sq with ti sufficiently 
close to So I and, moreover. 



p{x,ti) < -Ce{so - ti) 

in 

Brit^){zo)■ 


Now consider ip 

= 'ip{x,t) in such that 





o 

II 

on 

{\x- Zo 

= r{t)}, 


-A'iP_ =G{p_) 

in 

{ri < \x 

- zo\ < r{t)}, 

< 

'P+,t - = {ip+ - u)G(O) 

in 

{’’(i) < 

X- Zo\< r 2 }, 


Tp = p+{zo + r 2 V, t) + e{r 2 - r{t)) 

on 

{\x - Zo 

= r-2}, 


tp = p-{zo + riv, t) - s{r{t) - ri) 

on 

{\x- Zo 

= ri}. 


Note that (3.8) and the maximum principle for the elliptic equation in (1.12) yield that p{-,ti) < 
'tp{-,ti) in Br(ti){zo) \ Bri{zo)- 

The initial datum for ip in the parabolic phase at < = is given, for r = |a: — zo|, by 
{p{zo + r2V,ti) + e(r2 - r{t))){r2 - r(ti))"^(r - r(G)), 
so that, if r 2 is sufficiently close to r{ti), p{x,ti) < 'ip{x,ti) on Br^{zo) \ Bj.(ti)(- 2 o)- 
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Also note that, since ^ is a smooth strict supersolution of ( 1 - 12 ) with nonzero gradient at (xo,to)) if 
ri and r2 are sufficiently close to r, then and \Dtp-\ are close to (f) at (j/g, so) up to order e. 

Thus we can conclude that ■i/' is a supersolution of (1-12) in E if (ti,ri,r 2 ,e) is sufficiently close to 
(so,r, r, 0 ). 

Finally note that Um < </> < '0 on the parabolic boundary of E with (j) = tp at (yo, sg). 

The construction in the previous section yields a supersolution of (1-8) which converges uniformly 
to Ip as m ^ oo. By shifting, if needed, ^p down by a small amount, yields that Um crosses ip from 
below in the interior of E, a contradiction to the comparison principle of (1.1). 

□ 


Next we show that u and v coincide initiallly. For this we need (1.14) and (1.15). 

Lemma 4.2. Assume (1.14) and (1.15). Then, as t —)■ 0"^, 

(a) the sets d{u{-,t) > 0} and d{v{-,t) > 0} converge uniformly to Og in the Hausdorff distance, 
and 

(b) both u{-, t) and v{-, t) converge in fig to the harmonic function w in Qq with the zero boundary 
data, while outside o/f 2 g, they converge to the original initial data (uq)^ := p(l — pg). 


Proof. Assumptions (1.14) and (1.15) imply that, as m —> oo, there exist a sequence of positive 
numbers t 0 such that 


flo = {po = 1} c {u„,(-, 0) < 0} C {pg > 1 - Em), (4.4) 

and, in particular, the sets {Mm(',0) < 0}’s converge to ilg in the Hausdorff distance. 

We prove (a) first. Fix e > 0, a;g € Hg := {x : (i(x,Hg) < e} and let Hg(?/g) be the ball which touches 
dflo at Xg from outside of Hg. With a{t) := e/2 — Mt, we consider the radial function (p given by 


' -A0_ = G((/_) 

in 

{x : a{t) < |x - yg| < {G + 1)£} 

(/)_(•, t) = Mg 

on 

{|a:-2/o| = (G+ l)e}. 

(/(•,<) = 0 

on 

II 

O 

1 

4>t - vA(p ={(p- v)G{0) 

in 

{x : |x - Pg| < a(t)}, 

(/(•,0) = cg|x - ygp - Cg£2/4 

in 

the parabolic phase. 


with cg is chosen small enough so that </(•, 0 ) sits below «„(■, 0 ) in Bg/ 2 {yo); note that the support of 
</(-, 0 ) equals Bg/ 2 (yg). 

It follows from by Lemma A.l that \D(p+\{-,a{t)) is sufficiently large for some large M and for small 
enough t. Choosing such M and t we have \D(p_\ < on x = a{t), and, thus, (p a subsolution 

of (1.12) for a small time interval [0,tg]. Using the approximating solutions for the Um-equation 
constructed in the previous section as well as (4.4), we conclude that v > (p, and, hence, {v{-,t) < 0} 
is outside of Ba(t)(jJo) for 0 < t < tg. Since xg was arbitrarily and fig is compact, we find 

{v{-,t) < 0} C Hg for t < to = to{e). 


Letting £ —>• 0 gives 


Next we show that 


lim sup {v{-,t) < 0} C Hg. 
*->■ 0 + 


(4.5) 


flo C lim inf {u(-,t) < 0}. (4.6) 

*-> 0 + 

Fix B 2 r{xo) C Hg, let /(x,t) := t(|x —xgp —r^) and denote by p/ the pressure function corresponding 
to /. Then, as long as r is sufficiently small, 

ft - {mpf + iy)Af = (|x - xgp - r^) + 2nA{mpf + v)t\x\‘^ > -vG(pf) in B^(xg) x [0,1], 


Let 


^{x, t) . /(^T, ^)X|x—a:o|<r h{x^ 0X|tc—a:o|>r 5 
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where h{x,t) solves a heat equation in (M" — Br{xo)) x [0,1] with initial datum (uo)+ and zero 
boundary condition on dB^-ix^). Since (mo)+ = 0 in B 2 r{xo)^ it follows from the heat kernel estimate 
that, for small time. 


\Dh\{-,t) < < \Df\{-,t) on dBrixo). 


Thus g is a supersolution of (1.8) in K" x (0,to) for sufficiently small to depending on r, and, hence, 
for 0 <t<to, Br{xo) C < 0}. Since r is arbitrarily, (4.6) follows. 


In view of the facts that v < u and < 0} C {u(-,t) < 0}, the two inclusions above yield (a). 


Now we prove (b). Fix e > 0 and define 

rtf := {x : (i(a;,K" \ flo) > e} and ftg := {x : d{x,Qo) < e}. 

In view of (a), there exist 6 = S{e) > 0, to = to{e) > 0 and M such that for m > M and 0 < t < to 
the following holds: > S on dVlg, Um < —d in fly, and moreover 


—A/ = G(/) — e in ri/ and f = S on drif, 


and 

Let 


—Ag = G{g) + e in ilg and g = 6 on dVtg. 
(j){x,t) := —a{t)f{x) and il}{x,t) := —b{t)g{x), 


where 

a{t) := min[—1] and b{t) := max[(5“^e“”^^'^‘, Ij. 

Direct calculations then yield that, for sufficiently large m and any choice of to > 0, (j) and tp are 
respectively supersolution and subsoluton to the Um-equation inUf x (0,to] and rtg x (0,to]- Thus 
the comparison principle for the Um-equation and the choice of <5 and to yield 


■0 ^ '^m in X [Ojfo] and Um < (j) in flf x [0,to]- 

Letting m —> oo and using arbitrarily small e > 0, we conclude that the (um)-’s converge uniformly 
to the solution of the elliptic equation Dq with zero boundary data. 

Similar arguments apply to (um)+- Q 


To derive the convergence result for Um’s using the previous lemmata, we need to first show that, as 
m —)■ oo, the j/(l — Pm)’s stay nonnegative. 


Lemma 4.3. Assume (1.2), (1.3) and (1.5). Then 

\imsnppmix,t) < Pm in Qt, 

m—^oo 


and, hence, 


lim sup Pm A 1 in Qt ■ 

771—>-00 


Proof. Fix e > 0 and let C := max[pm(', 0),Pm +e]- Then G{pm +£) = — <5 < 0 for some <5 > 0. Thus 

f{x, t) := max[G — (m — l)G8t, pM + e] 

is a supersolution of (1.7) and Pm A / in Qt-, if Pm ^ f on the parabolic boundary of Qt- At t = 0 
this is guaranteed from the definition of G, and on the lateral boundary it follows from (1.3). Since 
/ = Pm + £ after t > 0(;^), we conclude that 

limsuppm(a:, t) < pm + e for all (x, t) S Qt- 

771—>-00 

The claim follows after sending e to zero. 
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□ 


We are now ready to show the main convergence result. 

Corollary 4.4. Let u and v be as in (4.1). Then: 

(a) u = V* is the unique viscosity solution w of (1.12) with the initial data uq. 

(b) As m ^ oo, the v{l — Pm) ’s and Pm’s converge uniformly respectively to b{w) and 1 — v~^b{w) 
in Qt- 


Proof. It follows from Theorem 4.1 and Lemma 4.2 that rt is a subsolution and t is a supersolution of 
(1.12) with initial data uq. Let w be the unique viscosity solution of (1.12) with initial data uq given 
in Theorem 2.5. Then u < w and w* < v, and, hence, 

u < w = (w*)* < V*. 

Since < it by definition, we have u = v* = w. A similar argument leads to t = w* = u*. 

Since, in view of Theorem 2.5 we know that b{w) is continuous, b{w) = 6(ic*) and thus 
b{u) = b{v) = b{w). 

Now we claim that 

b(v) < liminf*i/(l — pm) and limsup*p(l — pm) < b(u), 
which will yield the uniform convergence of the ;/(! — PmYs. 

To prove the first inequality in the claim above, we observe that Um < v{\ — pm) and, in view of 
Lemma 4.3, limsup^_^(^ pm < !> which gives liminfm->.oo v{l — pm) > 0. Then 

b(v) = max[0, liminf^Um] < liminf»p(l — pm)- 

To show the second inequality, suppose that, for some {xqjIq) G Qt, limsup*t/(l — Pm)ixo,to) > 0, 
since otherwise there is nothing to show. In this case, along a subsequence of m —)■ oo, we have 
liminf*pm(a;o; ^o) < cq < 1 and, thus, liminf*pm(aio, to) = 0. This yields that 

u{xo,to) = limsup*Wm(a;o,to) = limsup*p(l - pm)ixo,to), 
and the conclusion follows. □ 

The uniform convergence of Pm, or Corollary 4.4 seem to be the best one can do without further 
restriction on the initial data, even for the short time interval. In particular, when uq is a (strict) 
supersolution of the elliptic equation in (1.12), then the Um’s converge to the solution of (1.12) with 
the initial data strictly smaller than uq, that is one should expect, for sufficiently large m, almost 
discontinuous decrease of the pressure profile at the initial time for pm ■ 

5. Finite speed of propagation 

We derive some properties of the speed of propagation for both the negative and positive phases 
under certain assumptions. We know that the negative phase can nucleate in the interior of the 
positive phase, when u decreases from a positive value to become zero in finite time. We show here a 
quantitative version of this phenomena. Indeed when u becomes very small in an open neighborhood 
in the positive phase, then it turns into the negative phase in short amount of time. 

Lemma 5.1. Suppose u(-,to) < s in Br(xo) with 0 < e < r^. There exists a dimensional constant 
Cn such that 

C 

u{-,t) <Q in B^/^{xo) for ^ ^ Yo + + r^]- 

In particular, the interior o/{«*(•, t) < 0} coincides with {«*(•,<) < 0}. 
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Proof. Note that any solution to the parabolic equation 

Wt — vlS.w = {w — u)G(O) in Qx, 

which stays positive, serves as a supersolution to (1.12). This fact and the parabolic Harnack inequality 
yield 

u < Ce for a dimensional constant G > 0 in E := Br/2{xo) x [toi^o +1'^]. (5-1) 

Let f{x,t) := Ge + ^|a; — xqP — |uG(0)(t — to)- It is immediate that / is a supersolution to (1.12) 
and thus u < f. In particular, 


w < 0 in Br/iixo) X [ti,to + 


where 


ti 


to + 


4Ce 
uG(0) ■ 


□ 


The following two local estimates quantify the finite propagation property. First we show that the 
negative phase (the tumor) does not shrink too fast over time. Note that it can shrink at least 
temporarily, when the positive part of u has steep growth near the interface. In terms of the original 
density variable p this means that its is close to zero near the tumor boundary. 


Lemma 5.2. [Finite-speed Shrinkage] Suppose u < 1 in Qt- The negative phase {u{-,t) < 0} shrinks 
at most by tf^^ over a time period t. More precisely, there exits cq > 0 depending only on n, such 
that, whenever Br{xo) C {u{-,tQ) < 0} with 0 < r < cq and to > 0, then 

Br/ 2 {xo) C {u{-,t) < 0} for to<t < to 


Proof. Let h solve 

—Ah = G(h) in Bj.^ 2 ixo) with h = 0 on dBr/ 2 ixo). 

Then, since G(0) > 0, there exists G > 0 such that 

\Dh\ > Cr on dBr/ 2 {xo)- (5.2) 


Let Ip be the solution to 

( 'ipt- vAip = V'G'(O) in E := (M” - Br/ 2 {xo)) x [to, to + 

< 'tp = 0 on 5Br/2(xo) X [to,to + r^/2], 

[ = X{\-xo\>r} in iR^\Br/ 2 ixo)) X {0}. 

Using the properties of the heat kernel and scaling arguments we find, for some dimensional G > 0, 

\T>ip\i-,t) < G-e”’’ in dBr/ 2 {xo), (5.3) 

r 

and it follows that 

\D'4’\ < -e~^^'^ on the lateral boundary of E. 

Hence, if r is small, (p = (p^ — <p- defined by 

(p+{x,t) := ipix, t) in E, <p-{x, t) := h{x) in 5^/2(a;) x [to, to + 

is a supersolution of (1.12) in K" x [to, to + and the conclusion follows by comparing u with cp. 

□ 


Due to the possible nucleation of the negative phase, it is not feasible to estimate the speed of 
expansion of the negative phase without knowing the positive density distribution near a given point. 
The following lemma states that, if {xo,to) is on the interface and is nondegenerate near it, then 
the negative phase cannot expand too fast. 

Lemma 5.3. [Expansion bound] Suppose M(-,to) > 0 in B^rixo) with u{-,to) > r in B,.(xo). Then, 
for r < vq where tq depends only on n, we have u> Q in i?r/ 2 (a;o) x [to, to + 
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Proof. Without loss of generality we set xq = 0. Comparing u with a fundamental solution of the 
heat equation with initial data —X|a;|> 3 r eis a subsolution to (1.12), we find, for some C > 0 which 
again depends only on n, 

< Cor on dB 2 r{ 0 ) x [to,to +r^]. 

Next we construct a subsolution to (1.12) in i? 2 r x [to, to + by solving, for r(t) := r — y/t — to, the 
initial-boundary value problem 


(f-ix. 

II 

on 

-A(j). 

-{■,t) = G{cf) 

in 

dt(j)- 

> 

II 

1 

'0 

in 


= 0 

on 

?^+(’)^o) — SXBr-iO) 

in 


dB2r{0), 

B2r{0) - Br{t)i0), 
Uto<t<to+r3 {Br(t){0) X {t}) 
3S,.p)(0), 

{B2r\B.,(^t)) X {0}. 


The parabolic Harnack inequality for the heat equation gives \D(j)+\{-,t) > on 5^(4)(0). Thus 


■= W-i - (</*)- 


is a subsolution to (1-12) in B2r{0) x [to, to -|- r^], since, for a dimensional constant C, 

\Df>.\{r{t),t) < < C 2 ^f= < |i5<^+|(r(t),t). 

r y/t-to 

We can now conclude by comparing u with (j) in i?2r(0) x [to, to -I- r^]. 


□ 


Appendix A. Discussion on the proof of the comparison principle in [KP]. 

For fixed r > 0, we consider the open balls 

Br{x,t) := ^{y,s) :\x-yf+ \t-<r'^''^ and Bf{x):={y:\x-y\<r}, 
the space disk 

Dr{x) := Bf{x) X {0} = {{y, 0) : [a; - y| < r} , 

and the flattened set 

Er{x,t) := |(y,s) : \x - yf + \t - s\^ <r^|. 

Finally, we define the domain 'E.r{x,t) that is used in the definition of regularizations of solutions by 

Sj.(x, t) := Dr{x, t) + Er{x, t), 

where -I- is the Minkowski sum. Note that near its top portion at t = r, S(x,t) shrinks its spatial 
radius with the same rate. Finally, when (x, t) = (0, 0) we omit and we simply write 
We write next a simple version of the results stated in Lemma 3.13 and Lemma 3.16 of [KP] . 

Lemma A.l. (a) Let u > 0 solve Ut — vAu < 0 in a parabolic neighborhood of'E.r for some r € (0,1] 

and assume that u = 0 on S,.. There exists e > 0 and g{s) := —((s — t)+)^, where M is independent 

ofe, such that 

0 < u(x, t) < g{\x\) in {|x| < r -I- e} x [0,0 -I- rj. 

(b) Let V satisfy Vt — nAv > —1 in Sr(^, s) and v > 0 in Sr(^,s) for some (C, s) G K" x K. There 

exists f € C([0,r]), /(O) = 0, / > 0 on (0,r] and —>-00 as s ^ 0-1- such that 

v{x, t) > f{r — |x|) for {x : |x| < r} x [s -|- r/2, s -I- rj. 
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The first part of the lemma states that, if the support of u is S, then the spatial gradient of 
grows to infinity on the lateral boundary of S as t increases to r. The second part implies that, if the 
support of V is outside of S, then the spatial gradient of v{-,t) vanishes on the lateral boundary of S 
as t increases to r. 

Now we touch on the proof of the comparison principle. Since it parallels that of Theorem 3.1 in [KP], 
we only sketch it pointing out differences. 

We consider the regularized sub- and supersolutions 

Z{x,t) := sup u, W{x,t) := inf v, 

and argue by contradiction assuming that there is a finite first crossing time to defined by 

to := sup{r : for 0 < t < r}. 

It follows that there is a contact point at t = to between the free boundary of Z and W. 

Lemma A.2. We have: 

{Z{-,to) > 0} n {W(-, to) < 0} = d{Z{-,to) > 0} n d{W{-,to) > 0}. 

Proof. First note that, in view of Lemma 5.2, {M(-,to) > 0} can not expand discontinuously. The 
definition of Z then yields that, if xq lies in the interior of {Z(-,to) > 0}, then Z > 0 in Br{xo) x 
[to — t, to] for some sufhciently small r and t. It then follows from the definition of to and the strong 
maximum principle for the heat equation that {IF(',to) > 0} in the interior of {Z(-,tQ) > 0} and the 
claim follows. □ 

The next lemma corresponds to Lemma 3.19 in [KPj. Once it is shown, the rest of the proof is the 
same as in [KPj. To state it and sketch the proof, it is necessary to remind the reader the geometric 
properties of Z and W as in [KPj. By definition, for any yo G d{Z{-,to) > 0} there exists a “cylinder” 
Ez = Sj.( 2 :i,si) in {Z > 0} such that {yo,to) lies on the lateral boundary of Ez- On the other hand, 
for any j/o5{IF(',to) < 0} there exists an exterior “cylinder” E\y = Er{z 2 ,S 2 ) in {W < 0} such that 
(yoj^o) lies on the lateral boundary of Ew- 

Lemma A.3. Fix yo G d{Z{-,to) > 0} O d{W{-,to) > 0}, and let Ez and Eyy be as given above. 
Then 

to - r < si,S 2 < to + r. 

Proof. It follows from the definitions that {IF(-, to) ^ 0} has an interior ball of radius r which touches 
yo on its boundary. This fact and Lemma 5.2 yield that Si < to + Since Z <W prior to t = to, we 
also have S2 > to — r. If we show that S2 < to -|- r, then the ordering between Z and W implies that 
Si > to — r and, hence, the conclusion. 

Since si < to-|-r, E := S^njt < to} is nonempty and Z > 0 in E. Let B := S^njt = to}. The strong 
maximum principle as well as Hopf’s lemma imply that IF(-,to) > 0 in i? with nonzero gradient on 
the boundary of B. It then follows from Lemma 5.3 that the negative phase of W cannot expand too 
fast, that is S 2 < to -I- r. □ 
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